Current-voltage characteristic of narrow superconducting wires: bifurcation 

phenomena. 
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The current-voltage characteristics of long and narrow superconducting channels are investigated 
using the time-dependent Ginzburg-Landau equations for complex order parameter. We found out 
that the steps in the current voltage characteristic can be associated with bifurcations of either steady 
state or oscillatory solution. We revealed typical instabilities which induced the singularities in 
current-voltage characteristics, and analytically estimated period of oscillations and average voltage 
in the vicinity of the critical currents. Our results show that these bifurcations can substantially 
complicate dynamics of the order parameter and eventually lead to appearance of such phenomena 
as multistability and chaos. The discussed bifurcation phenomena sheds a light on some recent 
experimental findings. 
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The problem of the appearance of the resistive state 
in narrow superconducting channels when the current 
j exceeds the value of the critical current j c was the 
subject of intensive debates in the past[l|. The volt- 
age between two parts of the superconducting wire ap- 
pears when the superconducting order parameter (OP) 
is strongly suppressed near the phase slip center (PSC). 
The phase of the OP becomes very steep function across 
the PSC and finally demonstrates 2n jump exactly when 
the OP reaches o[l|, 0] ■ This process appears periodically 
in time leading to the oscillation of the OP and nonlin- 
ear current-voltage characteristic (CVC) with well pro- 
nounced steps of voltage when j reaches certain critical 
values Hi]. 

The PSC has been extensively investigated theoreti- 
cally and experimentally in the past (see Ref.[l[ and ref- 
erences therein). It is well established that if j > j% ~ 
jJ^/S, 2 the normal state is absolutely stable, here Ie is 
the penetration depth of the electric field into the super- 
conductor and £ is the superconducting coherence length 
Therefore if Ie » £, there is a wide interval of currents 
j c < j < ]2 where the superconducting state coexists 
with the normal state leading to the oscillating type of 
behaviour 0, 

The steps in the CVC are usually interpreted as points 
where new PSC penetrates to the wire [3, [3, 0|- It was also 
demonstrated that the time average of the voltage V and 
the frequency of oscillations oj satisfy the relation similar 
to the Josephson relation V cx u> The theoretical in- 
terpretation of the CVC is usually based on the statical 
model of the resistive stated, or a vei T similar dynam- 
ical theory of the resistive state [l], [lljj . All these models 
are based on the approximate time averaging of the time 
dependent Ginzburg-Landau equations (TDGLEs) over a 
period. Therefore the complex behaviour of the order 
parameter during one cycle is lost in these theories. An 
important step forward in the analysis of the TDGLEs 
has been made recently 13, HI- It has been found that 



the dynamic of the order parameter in the PSC is gov- 
erned by the ratio between the relaxation times of the 
TDGLEs [lEl El- It has been established that the CVC 
exhibits an S shape in the regime of the fixed voltage. 
Besides, it has been shown that in the limit of strong 
pair breaking effect due to interaction with phonons the 
CVC demonstrates strong hysteresis in the regime of con- 
stant current [l^. This behaviour appears because two 
different solutions may coexist at the same value of the 
current. All these results (see also 1J,[15[) demonstrate 
that different periodic and quasiperiodic in time solutions 
emerge with the changing of the voltage and current (l6| . 
However, some of the properties of these solutions are 
still not understood. The types of bifurcations when one 
solution is replaced by the other solution are not iden- 
tified. The theory still lacks of an accurate analytical 
description of the CVC and well defined and experimen- 
tally verifiable predictions. In this paper we present a 
detailed study of the CVC in the resistive state of super- 
conductors with different lengths. We reveal the types of 
bifurcations corresponding to the discontinuities on the 
CVC and derive an analytical expression for CVC when 
j — > j c . We predict the appearance of the frequency in 
the spectrum of the electromagnetic radiation generated 
by the current after the period doubling bifurcation. This 
frequency equals exactly to the half of the one before the 
bifurcation point. Besides, we identify the positions of 
the PSCs in the wire and prove that our predictions on 
the basis of TDGLEs are in excellent agreement with re- 
cent experimental results of Sivakov et al. [13] ■ 

Recent experimental studies have detected the spatial 
positions of PSCs in the different parts of the CVC[j.[l7l|. 
With the help of low temperature laser scanning mi- 
croscopy technique it was demonstrated that each quasi- 
linear part of the CVC has different number of the PSCs. 
When the current becomes larger then j c the first PSC 
appears in the middle of the wire. Further increase of the 
current leads to the second step in voltage. Two PSCs 
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appear symmetrically with respect to the center of the 
wire. Further increase of the current leads to the the ap- 
pearance of the third step in voltage and the third PSC 
is created in the center of the wire (see Fig. 2 of Ref . ) 

Another important effect was observed recently on the 
CVC of submicron BSCCO bridges^. When the tem- 
perature of the bridge is reduced bellow T c , characteristic 
steps on the CVC appear with clearly observed hystere- 
sis. The measurement of the voltage as a function of 
time at high temperatures and at the fixed current has 
demonstrated sharp switching between two metastable 
states. The voltage randomly jumps between two values 
(random telegraph noise) [181] . The average frequency of 
the switching is very strong function of temperature: de- 
creasing of temperature by 1 K leads to the decrease of 
the frequency of switching by several orders of magni- 
tude. Therefore the activation energy that corresponds 
to the switching is more then 10 4 K[18j, indicating the 
macroscopic nature of the barrier between metastable 
states. 

In our work, which involves numerical and analytical 
analysis of TDGLEs, we show that each step of the CVC 
corresponds to the bifurcation point where the solution 
looses stability and a new topologically different solution 
emerges. At j = j c the steady state solution looses stabil- 
ity as a result of the saddle- node homoclinic bifurcation 
19], which leads to the birth of a limit cycle of infinite 
period. With the further increase of j this limit cycle 
looses its stability leading to the period-doubling bifur- 
cation. This behaviour is universal and does not depend 
on the length L. The next bifurcation point is not univer- 
sal and depends on L. With the further increase of L we 
observe the effect of multistability. Two or more different 
limit cycles coexist at the same value of the current. 

The first TDGLE in dimensionless units has the form: 
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The second equation defines the stationary current 
through ID wire: 



j = -— + — (?/>* — - ip—) 
dx 2i dx dx 



(2) 



Here tp is the dimensionless complex OP, the distance 
is measured in units of the coherence length £ and 
time is measured in units of the phase relaxation time 
tq = 47r ^' T " , A is the penetration depth, a n is the nor- 
mal state conductivity, and c is the speed of light. The 
parameter u — ^ is a material dependent parameter, 
where is the relaxation time of the amplitude of the 
OP. The electrostatic potential <j> is measured in the units 
of </>o/27rcT0, where <fto = nhc/e is the flux quantum, e is 
the electronic charge and h is the Planck constant, the 
dimensionless current density j is defined in the units of 
0oc/87r 2 A 2 £. Nonequilibrium superconductors are char- 
acterized by an additional length scale Ie which describes 



the penetration of electric field. From Eqs. (|ll2[) it follows 
that Ie = ^/u 1 ^ 2 . The time-dependent Ginzburg-Landau 
equations for complex order parameter were derived for 
T < T c in late 60th d, [§]. The simplified version of 
the TDGLEs used in this work corresponds to the so- 
called gapless case with large concentration of magnetic 
impurities t s T c << 1, where r s is the relaxation time 
on magnetic impurities (To| . In that case parameter u 
is 12. More realistic generalized TDGLEs are derived 
in dirty limit Ti m T c « 1, where Ti m is impurity scat- 
tering time. These equations take into account inelas- 
tic scattering by phonons. In that case parameter u is 
5.79. The penetration depth of the electric field Ie is 
substantially enhanced if the inelastic scattering time 
is large enough. The enhancement of Ie may be ap- 
proximated in the case of gapless superconductivity if 
we assume that phenomenological parameter u is small 
u < 1 11]. This approximation is widely accepted in 
the literature and leads to qualitatively and often quan- 
titatively correct results above and below T c . Here we 
consider the superconducting wire of a length L with the 
following boundary conditions :p(— L/2) = p(L/2) = 1 
and d(j)(— L / 2) / dx — d(b(L/2)/dx — 0, where we define 
phase 9 and modulus p of the order parameter according 
to the relation ib — pexp(i9). The absence of the elec- 
tric field at the end of the wire determines the gradient 
of phase d6(-L/2)/dx = d9{L/2)/dx = j. 

We start with the analysis of the steady state solutions 
of the TDGLEs (HE]) which are defined as tp = (1 - 
k 2 ) 1 / 2 exp (ikx). The equation for k has the form j = 
(1 — k 2 )k, and for j < j c = 2/3^/3 it yields two roots 
fci < hi- The first solution corresponding to k\ is stable 
and the second (k 2 ) is unstable. 

When the current reaches its critical value j — j c these 
steady state solutions collide with each other (ki = k 2 ) 
and disappear, thus suggesting a saddle-node bifurcation. 
Indeed, the the stability analysis of the steady state solu- 
tion performed in Ref.[2| shows that only one Lyapunov 
exponent crosses when j = j c . Additionally, our nu- 
merical calculation revealed that at the moment of this 
bifurcation, a periodic solution (limit cycle) with very 
large period arises in the phase space of the system in the 
very close vicinity of the disappeared steady states. As j 
grows, the period of the limit cycle quickly decreases. All 
those facts evidence that at j — j c the systems undergoes 
the saddle-node homoclinic bifurcation jl9j ]. 

Note that the value of j c — 2/3\/3 is obtained for the 
wire of the infinite length. The boundary conditions for 
Eqs. (|ll2p defined above lead to the slightly different value 
of the critical current. Therefore the value of j c {L) should 
be determined from the solution of stationary equations 
(I1I2[) with the same boundary conditions. 

According to the theory of dynamical systems, the be- 
haviour of a system in the vicinity of the saddle-node 
homoclinic bifurcation can be described by the normal 
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form[l9| 

y = /3 + a((% 2 , (3) 

where y is a representative phase variable, /3 = 
-2 3 / 2 u(j - j c )/3(« + 2)i c , a(0) = -2 3 / 2 w(u + 2). To 
derive the normal form we have to follow the standard 
recipe Ref..l9]. It is convenient to rewrite the TDGLs 
(|H2p in the limit of L — > oo in terms of p and the gauge 
invariant scalar $ = <f> + 89/ dt and vector Q = 80/ dx 
potentials. 

udp/dt = <9V<9:r 2 +p(l-p 2 -Q 2 ) 
up 2 ® = d{p 2 Q)/dx 

j = -d<f>/dx + dQ/dt + p 2 Q (4) 

These equations have the steady state solution p = po = 
^/2/3, $ = 0, Q = 1/V3, j — jc, which becomes unstable 
and has one Lyapunov exponent Ao which crosses at 
j = j c . Expanding Eq. Q up to the second order in 
deviations from steady state solution, rescaling the time 
t/u t and projecting this equation on the direction 
of the eigenvector corresponding to Ao we obtain Eq. © . 
The normal form © allows to predict the period T of the 
periodic solution which corresponds to the limit cycle. 
Following to Eq.(32,23) of Ref.[20j we integrate equation 
© over y between ±yi where yi is of the order of 1. 
T = 2 tan" 1 (y ly /a(0)/l3)/y/a(0)p. Therefore, when (j- 
jc)/ 'jc "C 1 the period of oscillations is determined by the 
formula: 

To verify theoretical prediction we have performed ex- 
tensive numerical simulations of Eqs. © and ©, us- 
ing the fourth order Runge-Kutta method. The spatial 
derivatives are evaluated using a finite difference scheme 
of the second and the fourth order. The calculations were 
performed for the superconducting wire of three differ- 
ent length L = Lq, 2L and 4i (L /£ = 10.88) and for 
u = 1/2. The different length leads to a slightly different 
value of the critical current caused by the boundary con- 
ditions. The calculated critical currents are j c (L)/j c = 
1.016, 1.002 and 1.00019 respectively. Some properties 
of the solution are universal and do not depend on the 
length L. On the other hand, there are some features of 
the solution and the CVC which are not universal and 
depend on the length L. 

The first universal property of the solution is the 
dependence of the period of the oscillation on (j — 
3c{L))/j c {L) in the range of currents j c (L) <j< j c i(L), 
where j c i(L) is the first critical current where the limit 
cycle looses its stability. According to Eq.© T f=s 

9.62 [(j - j c (L))/j c (L)]~ 1/2 for u = 1/2. The results 
of calculations of the period for the oscillatory solution 
are plotted in Fig. 1 for all three lengths L together with 
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FIG. 1: The period of the solution as a function of current 
for three different lengths of the wire: L — Lo (triangles) , 
2Lo (circles), 4Lq (stars). The solid line represents the results 
of Eq.([S}. Arrows indicate the period-doubling bifurcation 
point jd(L). Insets represent projection of the limit cycle 
trajectory to the (p(0), E(0)) plane before (a) and after (b) 
bifurcation point j c i(L). Here E(0) is the electric field in the 
center of the wire. 

the analytical estimate, given by Eq.©. There is very 
good agreement of numerical results with Eq.© over 
more then four orders of magnitude in (j — j c (L))/j c (L). 
Since Eq.© describes the asymptotic behaviour of the 
period T in the limit (j — j c (L))/j c (L) — > the agree- 
ment becomes better when the current is closer to its 
critical value. 

The divergence of the period implies that the CVC 
near j c (L) should have nonlinear form V cx T" 1 oc 

[{j ~ Jc(L))/j c {L)] 1/2 (Fig.2). This type of dependence 
is due to the quantization rule, derived in[2l|, that con- 
nects the averaged electric field and the period of os- 
cillations. Indeed, as it is clearly seen from Fig.2, the 
voltage V shows the universal dependence in accordance 
with Eq.©. Note that similar CVCs have been reported 
earlier [12j. However, the authors did not recognize the 
type of bifurcation and did not derive the normal form 
for this bifurcation. As a result they were unable to get 
analytical expression for CVC in the vicinity of j c . An ac- 
curate solution of dynamical equations (|1I2I) and asymp- 
totic, predicted by Eq.©, do not confirm an approximate 
result V oc — l/ln((j —jc)/jc), derived in Ref.[l[. The 
square root behaviour of the CVC, however, is in agree- 
ment with the earlier result of Aslamazov and Larkin|22], 
who had considered the CVC of the short (L <§; £) su- 
perconducting bridge. 

For relatively small currents j c (L) < j < jci(L) the 
phase slip center appears periodically in time in the cen- 
ter of the wire (Fig.2) in the agreement with the experi- 
mental results [17]. This oscillating behaviour of the OP 
leads to the oscillating behaviour of the voltage and emis- 
sion of an electromagnetic radiation with the frequency 
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FIG. 2: The CVC of the channel with L = AL . Open 
circles represent the more complex limit cycles, which coexist 
with the main periodic solution (multistability) . Insets (a-d) 
show the space-time arrangements of the PSCs in different 
regions of the CVC. The inset (c) corresponds to open circles. 
Inset (e) shows the CVC for j < j c i(L) in comparison with 
analytical formula (see the text). 
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Further increase of the current leads to instability of 
the limit cycle. At j — j c i(L) the periodic solution be- 
comes unstable and a new limit cycle emerges as the re- 
sult of period-doubling bifurcation. This is clearly seen 
from Fig.l where a sudden increase of the period is ob- 
served at j = jci(L). It appears that this instability is 
universal and does not depend on the length L. 

The period doubling bifurcation could be easily recog- 
nized if we plot the cycle trajectories below and above 
j c i(L). It is important to plot trajectories only for the 
gauge invariant quantities to avoid unphysical effects. 
Therefore, in the insets to Fig.l we present the peri- 
odic solution in coordinates p[x — 0) and electric field 
E(x = 0) = —dtf>(x = Q)/dx. As clearly seen from this 
figure, the period-doubling bifurcation manifests itself in 
the transition from a single-loop (period- 1) limit cycle to 
a double- loop (period-2) limit cycle [23j |. 

This bifurcation also leads to the singular behaviour of 
the CVC. This singularity is more prominent for longer 
wires L — 4L where the jump of voltage is well pro- 
nounced, while for L = Lq it manifests itself only by the 
change of the slope (Fig. 3 a). The arrangement of the 
PSCs in space and time is plotted in Fig. 2. As it fol- 
lows from Fig. 2, at j — j c i(L) the position of one PSC is 
shifted up with respect to the center and the position of 
the next PSC is shifted down. Therefore, the period of a 
new limit cycle is twice of the period of the limit cycle be- 
fore the bifurcation point. Each period includes now two 
PSCs, one is slightly above and another is slightly below 
the center of the wire, which is well agreed with experi- 
mental observations 



the electromagnetic radiation generated by the current. 
It should be pointed out that the period doubling bifur- 
cation in the superconducting channel with fixed voltage 
was reported earlier [l^. Nevertheless there is substantial 



geometric difference between this two cases. In Ref.[16| 
the PSCs were always detected in the center of the wire, 
in contrast with our results. 

Note that the critical current for period-doubling bi- 
furcation jd(L) is not universal and strongly depends 
on the length of the channel L. According to Fig. 3 
tici(L) - j c )/j c = 0.125,0.012,0.0006 for L = L ,2L 
and 4Lo respectively. Indeed it demonstrates the expo- 
nential dependence on the length L, as it was proposed 
in Ref.[Hj]. 

The next bifurcation which is common for all stud- 
ied cases is the destruction of the periodic cycle. For 
all studied cases there exists the second critical current 
jc2{L)/j c = 1.44,1.027 and 1.001 for L = L a ,2L and 
4L respectively where the periodic solution looses its 
stability. This transition does not show any clear singu- 
larity on the CVC. In all studied cases instead of the peri- 
odic solution we have found oscillating, but non-periodic 
solution. The Fourier transform of the voltage contains 
broad features instead of (5-function like spectrum before 
the transition. The solution in the vicinity of j C 2{L) is 
characterized by long intervals of time where the solu- 
tion is almost periodic and topologically similar to the 
solution before the bifurcation point j < j C 2{L). Then 
the instability develops and solution moves far away from 
the periodic cycle. After few periods the solution is com- 
ing back again to a very similar oscillating trajectory 
and performs regular oscillations for many periods. To 
demonstrate this behaviour we calculate Poincare map 
for j C 2 < j < jc5 and L = 2Lq. The projection of 
the phase trajectory to the plane [p{— L/4), E(— L/A)) 
never crosses the space near the center of the trajectory. 
Therefore we chose Poincare section as a plane which 
crosses the center of this trajectory E(—L/4) = 0.04 for 
j/j c — 1.0278. In Fig. (4) we plot the projection of the 
Poincare map to p(0) and E(0) as a function of discrete 
time, when trajectory crosses the Poincare section. The 
long horizontal lines represent an almost periodic solu- 
tion. The time, during which solution stays near the 
periodic laminar orbit becomes larger and larger when 
7 — > irfiiL) + 0, as demonstrated in the inset of Fig. 4 



U Hg. According to Eq. (32,23) of Ref.dJ the time, 
which trajectory stays near the regular limit cycle is pro- 
portional to t ot (j ~ jc2(L))~ 1 / 2 . Inset in Fig. 4 demon- 
strates that t indeed increases as {j—j C 2(L))~ 1 / 2 indicat- 
ing that we deal with chaotic solution which is developing 
via the intermittence 



23] 



17]- As a result of period-doubling, 



a new frequency ll>2 — wi/2 appears in the spectrum of 



Further increase of the current for the case of L = 2L$ 
leads to an additional step of voltage (j C 3{L)/j c = 1.058) 
in the CVC. In that case two remote PSCs located sym- 
metrically with respect to the center of the channel are 
accompanied by the third one, which appears in the mid- 
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die of the wire (Fig. 2 d). This does not destroy chaos, 
although changes apparently its topological properties. 
The CVC in this region of currents have some non-regular 
behaviour (Fig. 3 b), which might be associated with ap- 
pearance and disappearance of the so-called "periodic 
windows" , areas of the parameter values, where limit cy- 
cles with a relatively large period become stable [23j ]. 

The next step in voltage at jci(L)/j c = 1.084 is asso- 
ciated with the creation of a new PSC, again without the 
change of the chaotic character of the solution. And fi- 
nally, the non-periodic orbit disappears and a limit cycle 
becomes stable in the relatively large window of currents 
Uc5{L)/j c = 1.098 < j/j c < 1.114). This bifurcation 
produces a voltage step on the CVC, and in contrast to 
the all previous cases, the voltage decreases at the bifur- 
cation point (Fig. 3 b). The corresponding limit cycle 
has four PSCs, each of which has different positions in 
the channel. Remarkably, such behaviour of PSC was 
previously observed in experiments [l7j ]. 

The situation for the longer channel L = ALo is slightly 
different. First, we have found that for j > j c i{L) there 
are few ranges of the current where other limit cycles 
coexist with the main periodic solution, similarly to the 
earlier report 12j. Therefore, at the same value of the 
current we observe different stable periodic solutions with 
their own attracting phase spaces. The multistability is 
found only in the case of long channel (L — 4Lq). In 
the case of L = Lq and 2Lq we were unable to find any 
manifestation of multistability. The region where chaos 
is developed via the intermittence for the long chan- 
nel L = ALq is much more narrow then in the case of 
L = 2Lq. The chaotic oscillations becomes very quickly 
unstable and new stable limit cycle with 3 PSCs appears 
a t j/jc > jcs(L)/jc = 1.0011 in contrast with the former 
case, where this solution has also an irregular chaotic be- 
haviour. This bifurcation manifests itself by the step in 
voltage on the CVC and is consistent with the observa- 
tion made in Ref. 17j . 

The effect of multistability revealed in our calculations 
may be related to the effects of switching between two 
metastable states, discussed in Ref. [181] . Indeed the exis- 
tence of two different limit cycles with the different topol- 
ogy of phase trajectories and different averaged voltage 
implies that relatively large thermodynamic fluctuations 
may switch from one trajectory to another. If the switch- 
ing is fast and the averaged time between switchings is 
large in comparison with the period of the limit cycle, 
the measured voltage will be averaged over many peri- 
ods. The switching requires overcoming the energy bar- 
rier between two limit cycles therefore the time between 
two switchings is exponential function of temperature. In 
the case of wide bridges the activation energy is macro- 
scopic and depends on the width of the sample. There- 
fore in the case when the time of measurements is larger 
than the period of the cycle and shorter than the average 
time between the switchings the fluctuations of the volt- 
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FIG. 3: The CVC of the wire a) L = L , b) L = 2L . Arrows 
indicate different bifurcation points. 



age looks like a telegraph noise and are similar to that 
observed in Ref. [181] . 

In conclusion, we have analyzed the current-voltage 
characteristics of the superconducting channel of a differ- 
ent length. We have proven that the steady state solution 
losses its stability as a result of the saddle-node homo- 
clinic bifurcation, which leads to a limit cycle of the di- 
verging period. This unambiguously leads to the singular 
V cx (j— jc) 1 ' 2 CVC. We have demonstrated that the sec- 
ond anomaly in the CVC is caused by the period doubling 
bifurcation. The main consequence of this bifurcation is 
the appearance of the new frequency in the spectrum of 
electromagnetic radiation which is equal exactly the half 
of frequency before bifurcation. The appearance of the 
third PSC is not universal. In short channels this bifur- 
cation appears when the solution has chaotic character. 
Increasing the channel's length L stabilizes a periodic 
solution with 3 PSC. In the intermediate range of cur- 
rents the periodic solution becomes unstable leading to 
the chaos, which is developed via the intermittence. The 
effect of multistability in the presence of thermodynamic 
fluctuations may be responsible for the switching between 
two periodic solutions and random jumps of voltage be- 
tween two averaged values. 
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